An NMR analog of the quantum disentanglement eraser 
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We report the implementation of a three-spin quantum 
disentanglement eraser on a liquid-state NMR quantum in- 
formation processor. A key feature of this experiment was 
its use of pulsed magnetic field gradients to mimic projective 
measurements. This ability is an important step towards the 
development of an experimentally controllable system which 
can simulate any quantum dynamics, both coherent and de- 
coherent. 
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One of the most intriguing effects in quantum mechan- 
ics is the "quantum eraser" Given an ensemble of 
identically prepared quantum systems, this effect is de- 
scribed by the loss or gain of interference in a subensem- 
ble that is determined by the outcome of the measure- 
ment of one or the other of a pair of noncommuting bi- 
nary observables, respectively. Thus a quantum eraser 
demonstrates the principle of complementarity without 
making use of the corresponding uncertainty relation. 
Quantum erasers have previously been demonstrated by 
optical [[| as well as atom (6| interferometry. In this 
Letter we use liquid-state NMR spectroscopy on pseudo- 
pure states 0] to demonstrate a novel "disentanglement" 
eraser, due to Garisto and Hardy ||, in which not only 
interference, but also entanglement, is lost or gained in 
a subensemble. An analogous quantum erasure proce- 
dure operating on a pair of Bell states has recently been 
used by Zeilinger's group to prepare an entangled three- 
photon state p. 

An important goal of our group is to design and build 
increasingly more powerful experimentally controllable 
devices capable of precisely simulating the dynamics of 
any quantum system with an equal or smaller Hilbert 
space dimension. Previously we have addressed the issues 
of coherent control (To| and pseudo-pure state prepara- 
tion [fill, and we are now developing methods for non- 
unitary quantum operations. The disentanglement eraser 
is of particular interest in this regard, because it allows 
us to show that NMR on pseudo-pure states is capable 
of reproducing even the decoherent dynamics associated 
with strong (projective) measurements on the members 
of the ensemble, which are needed to create or destroy 
entanglement in this eraser (cf. pT|]). 

It should be understood that the density matrices of 
the highly mixed macro-states involved in liquid-state 
NMR experiments can always be rationalized in terms of 
ensembles of disentangled micro-states . Consequen- 



tially, the ensemble-average observations on pseudo-pure 
states reported here do not prove the existence of the cor- 
responding entangled micro-states in the sample. Never- 
theless, because pseudo-pure states provide an equivalent 
representation of the underlying quantum dynamics, our 
experiments created exactly the same ensemble-average 
coherences that would have been observed if the same 
operations had been applied to the corresponding pure 
state ensemble, and this is sufficient for our purposes. 

In the disentanglement eraser two of the spins (qubits) 
in a GHZ (Greenberger-Horne-Zeilinger) state (l^Jl^l arc 
regarded as the components of a Bell state labeled by the 
state of an additional "ancilla" spin #1 (left-most), i.e. 



|Vohz> = ^(|0)|00) + |1)|11». 



(1) 



Assuming that the computational basis corresponds to 
the eigenvectors of the o~ z spin 1/2 operator, a projective 
measurement of the ancilla along z yields a mixture of 
separable states |00) and |11) labeled by the ancilla spin. 
This corresponds to the ensemble 
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where E± = ^(1 ± a\) expresses the density matrices 
|0)(0| (g> o-i ® oi = E\ and |1)(1| ® cri ® oi = E x _ in 
terms of the Pauli matrix a\ = a z ® a\ ® u\, where 
(7i = 1 0> <0 1 + 1 1 > < 1 1 and I = oi <g> oi <g> oi is the 8 x 8 
identity matrix. 

Alternatively, expressing the GHZ state in terms of the 
x eigenstates of the ancilla spin and the Bell states of the 
other two spins yields 



GHZ) 



MI<M), 



(3) 



where |ar±) = (|0) ± |1))/V2 and |<£±) = (|00) ± |11))/V2. 
Thus a projective measurement along the x-axis followed 
by a rotation of the ancilla back to z gives 



p x =E l + \cj> + )^ + \+E l _\^){cj>_\. 



(4) 



This is a mixture of complementary Bell states each la- 
beled by the state of the ancilla. Note that the partial 
trace over the ancilla in \i/jghz)('4>ghz\, Pz and p x are all 
equal to (|00)(00| + |ll)(ll|)/2 = E\E\+E 2 _E Z _, so that 
these states can be distinguished only if the information 
contained in the state of the ancilla is used. 

These effects were demonstrated by liquid-state NMR 
using as the qubits the three spin 1/2 carbons in a 13 C- 
labeled sample of alanine {C x O^ - C 2 H(C 3 H 3 ) - NH+) 
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in deuterated water. With decoupling of the protons fig ], 
this spin system exhibits a weakly coupled spectrum cor- 
responding to the Hamiltonian 

Hi n t = ^[v\<j\ + v 2 a\ + v 3 al 

+ \{Ji 2 a l y z + J2ao*<4 + Ji3<rl<T% (5) 

where the z/'s are Larmour frequencies and the J's the 
spin-spin coupling constants in Hertz. The experiments 
were carried out on a Bruker AVANCE-300 spectrometer 
in a field of roughly 7.2 Tesla, where the resonant fre- 
quency of the second carbon is 75.4713562 MHz. The 
frequency shifts of the other carbons with respect to the 
second are 9456.5 Hz for the first one and -2594.3 Hz for 
the third, while the coupling constants are J i2 — 53.7, 
J-23 = 34.6 and J13 = — 1.4 Hz. The T% relaxation times 
for the three spins are 21, 2.5 and 1.6 s, while the T 2 
times are 550, 420 and 800 ms, respectively. 

The pseudo-pure ground state was prepared from the 
thermal equilibrium state by the procedure summarized 
in Table 1, which uses magnetic field gradients (denoted 
by [V]) to dephase off-diagonal elements of the density 
matrix at strategic points along the way 0. Letting 
p eq = a\ + a\ + a z be the traceless part of the equi- 
librium density matrix (with all physical constants set to 
unity), the first two transformations in the table yield the 
state (V3/V32)al + (p\ + cr z )E^_ . Spins 1 and 3 may then 
be transformed into the state (V3/V32) (a\ + + &lal) 
by the efficient two-spin pseudo-pure state preparation 
procedure described in Ref. (llj (Eq. (47)), yielding the 
three-spin pseudo-pure ground state 



Pi 



= ^(^^-|) = (|000)(000|-I). (6) 



The logic network shown in Fig. 1 transforms this 
state into the pseudo-pure GHZ state, and then deco- 
here the ancilla as indicated. The GHZ state is ob- 
tained by rotating spin 2 (since J13 <§; J\ 2 , J23) to 
the x axis in the rotating frame with a 7r/2 y- rotation 
Ry(ir/2) = exp(— i(Ty7r/4), and then using it as the con- 
trol for a pair of c-NOT (controlled-NOT Jis)) gates to 
the other two spins. This pair of c-NOT's was imple- 
mented by the propagator 7V 13 I 2 = e lia *- a * )E - 7r/2 = 
e io- a .B_7r/2 e -io- a ,£!_7r/2 ( ensur i n g cancellation of the phases 
±i between the two factors). The overall sequence of 
transformations on the corresponding state vector is thus: 

|000)^5|0)(|0) + |1))|0)/V2 
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(|000) + |111))/V2= \ijj GHZ ) (7) 



Implementations of these operations in NMR by RF 
(radio-frequency) pulse sequences may be found in 
Refs. [[7|JTl Jliifl . The resulting pseudo-pure GHZ state is 
written in product operator notation as [|19| 



A V3 / 1J , ^2 3 , 13 

PGHZ - J^Pz^z + Vz^z + a z a z 



a x a y a y ~ a y a x a y) ) ( 8 ) 



^1^2 3 ^1^2 3 
~ a x a x a x - a y a y a x 



and has previously been studied by NMR in Refs. |20j|21j]. 

The coherences of pghz can be dephased, exactly as 
they would be by strong measurements of a z on all the 
individual systems in the ensemble, by means of magnetic 
field gradients similar to those used to prepare the initial 
pseudo-pure state (cf. jll]]). Specifically, a constant gra- 
dient V = dB z /dz applied for a period t along the static 
field axis z causes spin evolution under the Hamiltonian 
z l^\ Sj=i °z' where 7 is the gyromagnetic ratio of all 
the spins. This multiplies each coherence pki (fc ^ £) with 
a spatially dependent phase exp(— ijmke^zt/2), where 
rake is the coherence order p2[ (i.e. the difference in the 
z-component of the angular momentum in units of h be- 
tween the \k) and \£) states [^6)). Thus after such a gra- 
dient pulse the density matrix averaged over the sample 
volume satisfies pkt = for all k 7^ £ save for the zero 
quantum coherences (rrike = 0) . Because only one spin is 
dephased in the eraser experiments, only single quantum 
coherences are of consequence. 

This dephasing operation was made specific to those 
coherences involving transitions of the ancilla spin 1 by 
applying a ir pulse to the other two spins, after which a 
second gradient pulse of the same amplitude and dura- 
tion "refocuses" all the other coherences. At the same 
time it is necessary to also refocus the evolution under 
the internal Hamiltonian using n pulses selective for sin- 
gle spins. A sequence of RF and gradient pulses which 
accomplishes this is (in temporal order): 



Pi = M z -Ml-M,-Ml l 



[VI 
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(9) 



The corresponding effective (average) propagator is sim- 
ply e~ lVzo */ 2 . This dephases the ancilla spin in the same 
way as would a strong measurement of a\ on every mem- 
ber of the ensemble. To dephase the ancilla in the same 
way as would a strong measurement of u\ , one need only 
rotate the ancilla to the z-axis with a 7r/2 y-rotation 
R x _ v {-k/2), as follows: 



Pi 



[*/2]l. 



Pi 



(10) 



The ancilla is left along z for subsequent tomography. 
The results of P z and P x applied to pghz are 



Pghz 



4V2 



(ii) 



Pghz — ► ^\? Z (J Z + cr z a x a x - a z a y a y ). (12) 

These states were confirmed by full tomography (ItJ . Be- 
cause only the single quantum (mki — 1) coherences give 
rise to observable (dipolar) magnetization, it is necessary 
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to collect spectra not only following the dephasing oper- 
ation, but also following additional 7r/2 pulses selective 
for single spins, to rotate the uiki = and rriki > 1 co- 
herences, as well as the populations (diagonal elements), 
into observable single quantum coherences. Tomography 
was performed at the points of the procedure indicated 
in Fig. 1; the real parts of these four density matrices 
are shown in Fig. 2 (the imaginary parts were essentially 
zero) . 

The overall precision of quantum information trans- 
mission was quantified by an extension of Schumacher's 
fidelity ]25|], which takes into account not only system- 
atic errors, but also the net loss of magnetization due 
to random errors. This measure, called the attenuated 
correlation, is given by 



Tr(p the p exp ) 
Tr(p the p the ) ' 



(13) 



Here, p the is the measured pseudo-pure ground state 
Pini i transformed on a computer by the same sequence 
of unitary and non-unitary (measurement) operations 
to which it was subjected on the spectrometer to get 
p ex P . Note that, since Tr{p ex Pp ex P) < Tr{p the p the ), the 
Cauchy-Schwarz inequality implies that — 1 < c(p exp ) < 
1. 

The values of the correlation for each of the four to- 
mographic readouts were c(p**f) = 1 (by definition), 
<Pghz) = 0-88, c{pfP) = 0.92 and c(pfP) = 0.93. Al- 
though not included here for brevity, tomography on the 
state 1 0) ( 1 00) + |ll))/\/2 yields an attenuated correla- 
tion of 0.93, showing that spins 2 and 3 were entangled 
before the GHZ state was created. The increases in c 
are not unexpected, since the additional it and gradient 
pulses needed to mimic measurements on pghz are easily 
implemented with high precision, and the tomographic 
errors are estimated at ±5%. The leading candidates 
for the loss of correlation are pulse imperfections arising 
from RF field inhomogeneity, less than perfect RF pulse 
calibrations, and relaxation. The total time before data 
collection in the complete experiments was ca. 65 ms; the 
time required to prepare the GHZ state from the initial 
state was 21 ms. Since the T2 relaxation times of the 
spins varied from 400 — 800 ms, the net loss of magneti- 
zation due to relaxation in going from pj n j to pghz was 
3 — 5%. Thus, the additional loss due to pulse imper- 
fections etc. was about another 5% or so, confirming the 
high precision of the strongly modulated pulse sequences 
used in these experiments. 

In conclusion, we have used a three-spin liquid-state 
NMR quantum information processor to obtain a high- 
precision implementation of the dynamics, both coherent 
and decoherent, underlying Garisto and Hardy's "disen- 
tanglement eraser" , and have found that the experimen- 
tal results confirm the theoretically predicted conditional 
expectation values. This shows that we can judiciously 



and selectively render phase information macroscopically 
inaccessible in a way that precisely mimics the decoher- 
ence attendant on strong measurements. It should be 
noted that during this dephasing operation all interac- 
tions among the spins were refocused, and that only the 
macroscopically accessible information contained in the 
ancilla spin due to its earlier interactions with the other 
two was changed. This was nevertheless sufficient to con- 
vert the net triple-quantum coherence in pghz into a pair 
of double-quantum coherences, conditional on the state of 
the ancilla representing this information. Unlike previous 
eraser implementations, it was not necessary to explicitly 
read out this information in each member of the ensem- 
ble in order to see the conditional coherence, because 
this was done for us by the coupling of the ancilla to the 
other two spins while the spectra were being measured. 
This ability to convert macroscopic correlations into one 
another via well-defined microscopic (molecular) interac- 
tions is the essence of ensemble quantum computing. 
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FIG. 1. Logic network for the disentanglement eraser. Ini- 
tially, a pseudo-pure state on spins 1, 2 and 3 is created, 
pini = |000)(000| = E\E 2 + El. A f y-pulse is then applied to 
spin 2, followed by two controlled-not (c-NOT) gates to cre- 
ate the GHZ state (see text). Conditionality on the second 
spin being in the 1 1) state is represented in the network by a 
filled circle on its time line. Finally, the two complementary 
measurements, a\ and o\, are applied to spin 1. State tomog- 
raphy was performed to fully reconstruct the density matrices 
at the positions indicated. 
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4) [V]e i T5 (<T « +CT 2) e - i f CT >' e _i f ( a l+ a D 

TABLE I. This table shows the transformations used to 
obtain the initial state from the thermal state, p eq . 
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FIG. 2. Experimental density matrices reconstructed by 
tomography (in normalized units). The rows are enumerated 
in the standard computational basis, where for example 000 
represents the state label 000) . Although not shown, the 
columns are similarly labeled with the leftmost end repre- 
senting |000) and the rightmost end representing |111). Pini 
is the three-spin pseudo-pure ground state, and pghz is the 
pseudo-pure GHZ state. The last two plots are p z , which is 
pa hz after decohering spin 1 about the 2-axis, and p x , which 
is after decohering it about the a;-axis. (Note: pghz, pz and 
p x have been magnified by a factor of two for clarity). An 
amount of identity, chosen to optimize the input projection, 
was added to all experimentally measured density matrices. 
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